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Introducing polar coordinates around a point in Euclidean space reduces the 
Euclidian metric to the degenerate form 

(1) dr 2 +r 2 du 2 

where r is the distance from the point and duj 2 is the round metric on the sphere. 
If X is an arbitrary manifold with boundary, the class of conic metrics on X is 
modeled on this special case. Namely, a conic metric is a Riemannian metric on the 
interior of X such that for some choice of the defining function x of the boundary 
(x e C°°(X) with dX = {x = 0}, x > 0, dx ^Oon dX), the metric takes the form 

g = dx 2 + x 2 h on X° = X\dX, near dX. 

Here h is a smooth symmetric 2-cotensor on X such that ho = h\g x is a metric on 
dX. 

In fact a general conic metric can be reduced to a form even closer to (|l|) in 
terms of an appropriately chosen product decomposition of X near dX, that is, by 
choice of a smooth diffcomorphism 

[0, e) x x dX O c X, O an open neighborhood of dX . 

The normal variable in x € [0, e) is then a boundary defining function, at least 
locally near dX, and the slices F \ x =x have given diffeomorphisms to dX. Now 
such a product decomposition can be chosen so that 

(2) F*g = dx 2 + x 2 h Xl in x < e, 

where h x is a family of metrics on dX. 

This reduced form is closely related to the behavior of geodesies near the bound- 
ary. Up to orientation and parameterization there is a unique geodesic reaching 
the boundary at a given point p. In particular the normal fibration of X near dX 
given by the segments F([0, e) x {p}), p € dX, consists of geodesies which hit the 
boundary, each at the corresponding point p. 

We shall discuss here the behavior of solutions to the wave equation 

(3) (D 2 - A)u = on R x X° 

when X is endowed with a conic metric, A is the associated (positive) Laplacian 
on functions, and D t = —id/dt. For simplicity we take X to be compact. It is only 
really important that dX be compact. 

Our primary concern is to describe the phenomenon of the propagation of singu- 
larities for solutions to (||). To do so it is necessary to understand the behavior of 
solutions in a way related to the functional analytic domain of A. For the moment 
we simply say that we are dealing with 'admissible' solutions. This condition is 
explained further below. 
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In the interior of X the propagation of singularities, described precisely in terms 
of the notion of wavefront set, was treated in detail by Hormander (0). We para- 
phrase Hormander's result here as 

"Singularities travel along null bicharacteristics, which in the case 
of the wave equation project to time-parameterized geodesies." 

Thus, in the microlocal sense of singularities described by the wave front set, a 
bicharacteristic segment, which covers a light ray, either consists completely of 
singularities for a given solution or the solution has no singularity along it. 

This quite adequately describes the propagation of singularities except where a 
light ray hits the boundary at some point and at some time. Here a 'splitting' of 
singularities will usually take place. This is generally called a diffractive effect. The 
contrapositive of this effect can be succinctly stated as follows: 

"If no singularity reaches the boundary at time i then no singular- 
ity leaves at time t." 

The point here is that the regularity along any one of the 'radial' rays leaving 
the boundary at a given time is related, in general, to the singularities on all the 
incoming rays (although there are two separate components, as described below) 
arriving at the boundary at that time. Thus, even if singularities arrive at the 
boundary at time t along just one ray, they will in general depart along all rays 
leaving the boundary at time t. 

There are, however, some important exceptions to this general spreading of sin- 
gularities. For instance let X be a conic manifold with 'trivial' conic metric defined 
by the blowup of a point in a smooth Riemannian manifold. In this case, because of 
Hormander's theorem on interior singularities, the singularities are carried outward 
only on the one ray continuing the incoming ray in the original manifold. 

For a general conic metric there is a similar notion of the 'geometric continuation' 
of an incoming geodesic which hits the boundary. For a trivial conic metric obtained 
from a blowup, the boundary metric h^ is the standard metric on the sphere. The 
geometrically related incoming and outgoing rays hit this sphere at antipodal points; 
these can also be thought of as the points separated by geodesies of length 7r on 
the unit sphere. In the case of a general conic metric we mimic this by defining the 
relation 

(4) G(p) = 

{q £ dX; 3 a geodesic in dX for ho of length ir with end points p, q}. 

In general of course, G(p) is not smooth, but generically it is a hypersurface with 
Lagrangian singularities; it is always the projection of a smooth Lagrangian relation. 

A geometric refinement of the diffraction result is obtained by considering the 
order of singularity with respect to Sobolev spaces and an additional 'second mi- 
crolocal' regularity condition. For simplicity suppose that the (admissible) solution 
u is singular only near dX and only near a single incoming ray hitting the bound- 
ary at time i and at the point p. In the past (for t < t) we may suppose that the 
solution is locally in some Sobolev space H s . Suppose further that the singularities 
of the solution are not too strongly focused on dX insofar as tangential smoothing 
raises the overall regularity, that is, for some k,£ > 0, 



(5) 



(A + iy k u e Hgf int<t near dX. 
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Figure 1 . Regularity occurs at outgoing rays from points not tt- 
related to p. 

Under these two assumptions and the additional requirement that 

(6) 0<£< 2' 

we obtain the following 'geometric theorem'. 

"If an admissible solution is singular only near an incoming ray 
arriving at dX at time t and (||) and (Q) hold, then on outgoing 
rays with initial point in the complement of G(p), 

(7) u G H^^ e V e > in t > i near dX. 

When slightly generalized, as described below, this result applies to the funda- 
mental solution 

sin tVA 

with pole close to dX and with i < ^y^- The diffractive theorem merely tells us 
that if the pole is specified at (x,p) at t = 0, then singularities cannot emanate 
from dX except at time t — x. On the other hand, while 'strong' singularities can 
emanate from all points in G(p), the geometric theorem tells us that the solution is 
microlocally more regular on rays starting from dX at t = x but with initial point 
outside G{p). 

In the special case in which the metric g takes precisely the 'product' form 

(8) g = dx 2 +x 2 h(y,dy), 

near the boundary, Cheeger and Taylor [Q, || have given an explicit analysis of 
the fundamental solution constructed by separation of variables. (See also the 
discussion by Kalka-Menikoff Q.) They show a stronger form of the result above, 
including a more precise regularity estimate. A version of the results of Cheeger- 
Taylor has been established in the analytic category by Rouleux |l3| . Lebeau |Q || 
has also obtained a diffractive theorem in the setting of manifolds with corners in 
the analytic category. 

The first author acknowledges partial support from NSF grant DMS9625714. 
The second author was partially supported by an NSF VIGRE instructorship at 
Columbia University and by an NSF postdoctoral fellowship, and is grateful to 
Alberto Parmeggiani and Cesare Parenti for helpful conversations. 
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1. FRIEDRICHS EXTENSION 

To describe the admissibility condition, near the boundary, for solutions to (|^) 
we first describe the domain of the Laplacian for a conic metric. We take the 
Friedrichs extension of A. By definition, A is associated to the Dirichlet form 

(9) F(u,v)= [ (du,dv) g dg , u,v eC™(X°) 

Hence, dg is the metric volume form; in this case 

dg = Lpx n ~ 1 dx dh near dX, n — dimX , (p S C°° , <p > . 

The inner product in (||) is that induced, by duality, by the metric on T*X°. 
Following Friedrichs we define, 

DiA 1 ' 2 ) = clos {C c °°pn w.r.t. F(u,u) + Nil,*} , 

whenever X is a compact conic manifold with boundary of dimension n > 2. 

This is a Hilbert space with dense injection D{A 1 / 2 ) «-» L 2 (X) so there is a dual 
injection L 2 (x) <-> (D(A^ 2 ))'. The natural operator A : D{A 1 / 2 ) -> (^(Zi 1 / 2 ))' is 
determined by 

(2\u, ^) L? = F(u, ^)V«,^ DiA 1 ' 2 ). 

Then the Friedrichs extension of A is the unbounded operator with domain 

D(A) =|u£ D(A 1/2 ), Au G . 

In this case, it is a self-adjoint, non-negative operator with discrete spectrum of 
finite multiplicity. This allows its complex powers to be defined by reference to 
an eigenbasis. The real powers are isomorphisms off the null space, which consists 
precisely of the constants. Each of the powers is therefore a Fredholm map 

A s : D{A S ) -> L 2 g {X) VseR. 

with null space the constants and range the ortho complement of the constants. The 
domains form a scale of Hilbert spaces, and 

D(A S ) ^ D(A t ) is dense V s > t 

with D(A°) =L 2 g {X). 

For our purposes it is also important to note when the domains consist of ex- 
tendible distributions, i.e. those dual to C°°(X). This is the case only for s > — *f 
and more precisely 

(10) C°°(X) ^ D(A S ) ^C-°°{X) 

are dense inclusions ior —j < s < j. The limits of this range correspond to the 
occurrence of formal solutions of Au = 0. 

2. Wave group 
The Cauchy problem for the wave equation 

(Df - A)u = 0, onlxl" 
u\ t =o = u , D t u\ t =o = ui 
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has a unique solution 

ue C°(R; J D(Z\ 1 / 2 ))nC 1 (R;L2( x )) 
V (u Q , Ul )eS = Sx = D(A^ 2 ) ffi L%X) . 

These 'finite energy solutions' are the main object of study here. More generally, 
with the equation interpreted in C _00 (K; D(Ai the Cauchy problem has a 
unique solution 

u G C°(M; D{Ai)) n C X (K; 

V(uo,ui)ef a = D(^*)©i?(2l4-4). 

The regularity hypothesis on the solution can be weakened to 

u G if oc (R; D(A? )) n P^K; I>(4*-*)) 

without changing the unique solvability. 

Notice that these calculations are consistent under decrease of s. Furthermore, 
partial hypoellipticity in t shows that the solution to ( pi] ) satisfies 

(13) u G i? l ; c fe (M;L>(Z\5+l)) Vfcet. 

An admissible solution to the wave equation is one that satisfies 

(14) u G P£ C (R;P(4*)) for some q,p G R 

with @ holding in J C /" 1 ~ P_2 (R; D(A% ~ 1 )). Such a solution automatically satisfies 
( fl3| ) for some s. 

These statements can be reinterpreted in terms of the wave group 

3. HORMANDER'S THEOREM 

Let M be a manifold without boundary. The wave front set of a distribution 
u G C~°°(M) is a closed subset of the cosphere bundle 

WF (it) c S*M. 

It may be defined by decay properties of the localized Fourier transform, or the 
FBI (Fourier-Bros-Iagonitzer) transform, or by testing with pseudodifferential op- 
erators. The projection 7r(WF(u)) C M is exactly the C°° singular support, the 
complement of the largest open subset of M to which u restricts to be C°°. 

A refined notion of wavefront set is the Sobolev-based wavefront set, denoted 
WF S ; this is a closed subset of S*M, where now the projection is the complement 
of the largest open subset of M to which u restricts to be H s . 

If it satisfies a linear differential equation, Pu — 0, then 

WF(tt) c E(P) C S*M 

when S(P) is the characteristic variety of P, the set on which its (homogeneous) 
principal symbol, p, vanishes. 

If p is real then the symplectic structure on T*M, or the contact structure on 
S*M, defines a 'bicharacteristic' direction field Vp on S*M, tangent to S(P). The 
integral curves of Vp are called bicharacteristics; those lying in £(P) are called null 
bicharacteristics. 
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Theorem 1 (Hormander). Let P be a (pseudo)- differential operator with real prin- 
cipal symbol. If Pu = then WF(u) C £(P) is a union of maximally extended null 
bicharacteristics. 

The same result also holds with WF replace by WF S for any s. 

In our case, M = Rx X° so T*M = T*R x T*X°. The principal symbol of the 
d'Alembertian is t 2 — | • \ 2 gl where r is the dual variable to t and | • \ g is the (dual) 
metric on T*X°. Then 

£(P) = £+(P) U S_(P) C S*M 

when £r(P) = R x S , *X° arc the disjoint parts of S(P) in r > and r < 0. In 
this representation of S(P) the null bicharacteristics are geodesies on X°, lifted 
canonically to S , *X°, with t as affine parameter. Thus, for the wave equation 
over X°, Hormander's theorem does indeed reduce to the informal propagation 
statement desribed above. 

Combined with standard results relating the singularities of the solution to sin- 
gularities of the initial data, Hormander's theorem applied to the wave equation 
on a conic manifold to yields complete information on the behavior of singularities 
except along bicharacteristics lying above geodesies which hit the boundary. 

4. DlFFRACTIVE THEOREM 

On parametrized geodesic segments with an end point on the boundary, the defin- 
ing function x is either strictly increasing or strictly decreasing near the boundary. 
For each sign of t and for each t £ R the bicharacteristics covering such geodesies 
which hit the boundary at t = t and along which t is increasing (resp. decreasing) 
as x decreases, form a smooth submanifold of {t < t} x X (resp. {t > t} x X). We 
denote these 'radial' surfaces (near dX) by 

R±j{t) and R±,o(i) C S(P) 

where ± is the sign of r and /, O refers to whether these are 'incoming' or 'outgoing' 
and hence, equivalently, whether they lie in t < t or t > i. 

Theorem 2 (Diffractive theorem). If u is an admissible solution to (|3|) then for 
any i e R, s e R, o = ±, 

R aJ (t) n WF s (u) = R^ (t) n WF s (u) = 0. 

Here, WF S («) is the wave front set computed relative to the Sobolev space H s , 
locally in the interior. 

This is a precise form of the diffractive result described informally above. Notice 
that the singularities for different signs of r are completely decoupled. This does 
not, however, represent any refinement in terms of propagation along the underlying 
geometric rays, since all geodesies are covered by bicharacteristics with r fixed and 
of either sign. 

The proof of this result is discussed briefly below in §]7j. 

5. Geometric theorem 

Consider a geodesic on X which hits the boundary at a point p S dX. An open 
set of perturbations of the geodesic, meaning geodesies starting near some interior 
point on the geodesic and with initial tangent close to the tangent to the geodesic, 
will miss the boundary. A limit of such curves as the perturbation vanishes consists 
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Figure 2. A sequence of geodesies nearly missing the boundary, 
and the three segments to which they limit. 

of three segments. The first is the incoming geodesic segment. The second is a 
geodesic segment in the boundary, of length ir. The third is the outgoing geodesic 
from the end point of the boundary segment, which is therefore a point in G(p) as 
defined in (^) (see Figure ||). Thus it is reasonable to suppose that, amongst the 
outgoing bicharacteristics leaving the boundary at time i, those with initial points 
in G(p) will be more closely related to an incoming bicharacteristic with end point 
p arriving at time i. We call these the geometrically-related bicharacteristics (or 
geodesies). 

For instance, if there are incoming singularities on a single ray the singularities 
on the 'non-geometrically-related' outgoing bicharacteristics might be expected to 
be weaker than the incoming singularity. However, this is not in general the case. 
To obtain such a geometric refinement of the diffraction result we need to impose 
an extra 'nonfocusing' assumption. 

Theorem 3 (Geometric theorem). Let u be an admissible solution to (|^) and let 
a = ±. Suppose that R a j(t) n WF S u ~ near dX . Suppose additionally that for 
some k and < £ < ^ 

(16) WF s+e (l + A ax )- k unR a j(t) = Q). 

For any < r < £—1/2, if no incoming bicharacteristic hitting the boundary at time 
i at a point in G(p) with sgnr = a is in WF s+r u, then the outgoing bicharacteristic 
with initial point p € dX and sgn t = a is not in WF s+r u either. 

Thus the additional assumption ( |l6| ) allows regularity on outgoing rays to be 
deduced from regularity in the incoming geometrically-related rays up to the cor- 
responding level above 'background' regularity. 

As already noted, this result may be applied to the fundamental solution with 
initial point near the boundary. If the initial pole of the fundamental solution is 
sufficiently close to the boundary then there is a unique short geodesic segment 
from it to the boundary, arriving at a point p. If i is the length of the segment then, 
provided t is small enough, ( |l6| ) holds with s < — ^ + 1 for any I < ^j^. It follows 
that on R±,o{i), the outgoing set, the fundamental solution is in H~ e , for all e > 0, 
microlocally near the non-geometrically related rays, those with end point not in 
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G(p), whereas the general regularity is iJ~Tr +1 ~ e for all e > 0. This is a gain of 
'nearly' \ — 1 derivatives over the background regularity. 

In this way we extend part of the result of Cheeger and Taylor jj], |^] in the 
product case (||) to the general conic case. Note, however, that inspection of the 
fundamental solution constructed in Q] reveals the diffracted wave to be (n— 1)/2— e 
derivatives smoother than the incident wave; here, we only obtain (n — 2)/2 — e 
derivatives of improvement, hence our results are not sharp in the product case. 
Even in the non-product case, we conjecture that the better estimate should hold. 

6. Spherical conormal waves 

Around a given point q in a compact Riemann manifold there are 'spherical' 
conormal waves which are singular only on the spherical surfaces r = ±t, for small 
t of both signs. These just correspond to conormal data at t = at the (Active) 
cone point q. An important example is the fundamental solution, in which case 
the result follows from Hadamard's construction. In the more general case of a 
conic manifold with boundary there are similar contracting, and then expanding, 
conormal waves. 

Theorem 4. // u is an admissible solution near dX and t — which is conormal 
to t = —x for t < then it is conormal to t = x, near the boundary, for small 
t > 0. 

These conormal solutions to the wave equation in the general conic case are at 
the opposite extreme to those considered in the Geometric Theorem above. Namely, 
they are already smooth in the tangential variables, so no tangential smoothing in 
the sense of (||) is possible. Further analysis of the structure of these waves shows 
that the principal symbols undergo a transition at x = 0, the boundary, given by 
the scattering matrix for the model cone with the same boundary metric. Since this 
scattering matrix should have full support in general, this provides counterexam- 
ples to any extension of the geometric theorem in which the tangential smoothing 
condition is dropped. 

7. Methods 

The basic method we use is microlocal, but non-constructive. It is a direct 
extension of one of the proofs by Hormander of the interior propagation theorem. 
This 'positive' commutator method is itself a microlocalization of the energy method 
for hyperbolic equations. In it a 'test' pseudodifferential operator, A, is applied to 
the equation and the essential positivity of the symbol of the commutator i [P, A] 
gives a local regularity estimate on the solution. 

To extend this method to cover behavior of solutions near the boundary we 
replace the ordinary notions of wavefront set, pseudodifferential operators and mi- 
crolocalization with versions appropriately adapted to the geometry. When con- 
sidering the Laplacian itself on the manifold with boundary with conic metric, the 
appropriate notion is that of a weighted &-pseudodifferential operator (see Jl2|). 
This for instance allows the precise description of the domains of the powers of A 
which is used at various points in the argument. 

However, for the wave operators for the conic Laplacian the appropriate notion 
corresponds to the 'edge' calculus of pseudodifferential operators discussed origi- 
nally by Mazzeo || , arising from a filtration of the boundary (see also Schulze Q ) . 
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In this case, the manifold with boundary is X x R and the fibers of the boundary 
are the surfaces t = const. Thus t is the base variable of the fibration. 

To the edge calculus of pseudodifferential operators, given by microlocalization 
from the differential operators generated by xD x , D y (where the j/'s are tangential 
variables) and xD t , we associate a notion of wavefront set. We can prove the 
propagation theorem analogous to that of Hormander for this 'edge' wavefront set. 
However, in this new sense, D\ — A is not globally of principal type but rather has 
two radial surfaces. These correspond to the end points of bicharacteristics arriving 
at, and leaving from, the boundary. At these surfaces there are restrictions on the 
propagation results, very closely related to those for scattering Laplacians in p0| . 

These propagation estimates form the basis of both the diffractive and geometric 
theorems. In the former we combine the estimates with a variant of the one- 
dimensional FBI transform, scaled with respect to the normal variable x. This 
reduces the diffractive result to an iterative application of a uniqueness theorem for 
the Laplacian on the model, non-compact cone. 

To obtain the geometric theorem, showing that the outgoing singularities on 
non-geometrically related rays are weaker than the incoming ones, we use a division 
theorem. The additional hypothesis of microfocal tangential smoothing is shown to 
imply that the solution actually lies in a weighted Sobolev space with a higher x 
weight (hence more 'divisible' by x) than is given, a priori, by energy conservation. 
This allows the microfocal propagation results indicated above to be pushed further 
at the outgoing radial surface and so yields the extra regularity. 



8. Applications and extension 

The propagation of singularities results of the type discussed above should allow 
estimates of the spectral counting function as shown originally by Ivrii ( Q , see also 
llj and |). 

We expect these methods to extend to more complicated geometries, including 
manifolds with corners and iterated conic spaces. 
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